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In a nuclear medium, spin-dependent forces cause the nucleon spins to fluctuate with a rate r CT . We 
have previously shown that as a consequence the effective axial-current neutrino-nucleon scattering 
cross section is reduced. Here, we calculate this reduction explicitly in the perturbative limit r CT <C T. 
By virtue of an exact sum rule of the spin-density structure function, we express the modified 
cross section in terms of the matrix element for neutrino-nucleon scattering in the presence of 
a spin-dependent nuclear potential. This representation allows for a direct comparison with and 
confirmation of Sawyer's related perturbative result. In a supernova core with a typical temperature 
T = 10 MeV, the perturbative limit is relevant for densities p < 10 13 g cm -3 and thus applies around 
the neutrino sphere. There, the cross-section reduction is of order a few percent and thus not large; 
however, a new mode of energy transfer between neutrinos and nucleons is enabled which may be 
important for neutrino spectra formation. We derive an analytic perturbative expression for the 
rate of energy transfer. 

PACS numbers: 97.60.Bw, 13.15.+g, 14.60.Lm, 95.30.Cq 



I. INTRODUCTION 

Neutrino scattering rates in a medium differ from those 
taking place in vacuum. It is well known that spatial 
correlations between the locations or spins of the target 
particles can reduce or enhance the average effective scat- 
tering cross section. For example, the anticorrelations 
caused by the Pauli exclusion principle are straigthfor- 
ward to include. Even in a nondegenerate medium, cor- 
relations are induced by forces between the targets such 
as the Coulomb force which thereby causes electromag- 
netic screening effects jjj . Similarly, in a nuclear medium 
the spin-dependent nature of the nucleon-nucleon inter- 
action may cause nonnegligible "pairings" of the nucleon 
spins and thus a reduction of the axial-current neutrino- 
nucleon scattering rate ||. 

We presently study a less familiar cross-section modi- 
fication which is caused by temporal fluctuations rather 
than spatial correlations. The main idea is that the 



neutrino scattering process takes a certain amount of 
time. If the energy transfer is u, the weak probe can- 
not "resolve" those temporal changes of the target con- 
figuration which take place on a time scale faster than 
about 1/oj. For example, the target nucleon spin may flip 
"during" the neutrino-nucleon collision and thus "cancel 
itself." In linear-response theory, this effect is formally 
described by the frequency dependence of the nucleon 
dynamical spin-density structure function, which in the 
relevant limit amounts to the Fourier transform of the 
autocorrelation function of a single nucleon spin. In the 
absence of interactions the nuclear spin and thus its au- 
tocorrelation function is constant. In the presence of a 
spin-dependent random force the initial spin direction is 
forgotten, causing the spin autocorrelation function to 
decrease to zero for large times. Loosely speaking, then, 
for small u (large "duration" of the collision) the weak 
probe sees a reduced average target spin and thus scatters 
less efficiently. 
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A complete treatment should simultaneously include 
spin fluctuations and spin-spin correlations, and presum- 
ably spin waves as well. All these effects are embedded 
in the dynamical spin-density structure function, which 
in general has multiple isospin components. However, in 
contrast to spin-spin correlations, spin fluctuations occur 
even when there is only one nucleon — provided that its 
spin is jiggled around. This is a multiple-scattering effect, 
not a many-body phenomenon. In certain circumstances 
a pion condensate H or the walls in the nuclear bubble 
phase Q may be the dominant cause for nucleon spin 
fluctuations so that it is not entirely academic to study 
spin fluctuations independently from spin-spin correla- 
tions. 

Collision-induced changes of particle velocities or spins 
cause the bremsstrahlung emission of photons, neutrino 
pairs, or axions. According to the Landau-Pomeranchuk- 
Migdal (LPM) effect HQ the low-energy part of the ra- 
diation spectrum is suppressed if multiple interactions 
destroy the temporal coherence of the source. The spin- 
fluctuation effects studied here are analogous, except that 
it is the neutrino scattering rate that is being reduced. 
While the LPM effect is usually discussed for vector- 
current processes and thus for velocity fluctuations, in 
the case of axial-current processes in nonrelativistic nu- 
clear matter the spin fluctuations are more significant. 
We note that temporal fluctuations do not occur for a 
conserved quantity such as the charge of a particle. The 
vectorial nucleon quantity that does fluctuate due to col- 
lisions is the velocity, which in the nonrelativistic limit 
is small. Therefore, in this limit multiple-scattering ef- 
fects are not important for vector-current neutrino inter- 
actions Still, because in vacuum the nonrelativistic 
neutral-current neutrino-nucleon scattering cross section 
is a = (Cy + 3C^)G|£^/47r, any modification of the 
axial-current part strongly affects the total rate. 

The importance of multiple scattering is quantified by 
the spin fluctuation rate T a which roughly represents 
the inverse of the time required for the nucleon to for- 
get its initial spin orientation. This effect is important 
if r CT is of order the typical energy of the weakly inter- 
acting particles which scatter off, or are emitted from, 
the medium [a], i.e. for r CT > T. One can easily es- 
timate (Eq. Il5l below) that in a supernova (SN) core 
with a temperature of order lOMeV this "high-density 
case" obtains for p > 10 13 gcm~ 3 . Because densities 
as large as 10 15 gcm~ 3 are encountered in a SN core, 
quantities like the neutrino opacity or the axion emis- 
sivity are impossible to calculate in a purely perturba- 
tive way which is based on the assumption that aver- 
age scattering or emission rates are the incoherent sum 
of single-scattering events. Interaction rates calculated 
in the "vacuum limit" are fundamentally flawed for the 
conditions of a SN core. 

To extract meaningful estimates for weak interaction 
rates one must take recourse to the more general princi- 
ples of linear-response theory. With our collaborators we 
have begun to develop this perspective in a series of pa- 



pers |^,p|-pH| . We have argued that the neutrino opacities 
or axion emissivities can be estimated by virtue of a phe- 
nomenological ansatz for the spin-density structure func- 
tion which incorporates certain limiting cases, notably 
the low-density one, and which satisfies certain general 
principles, in particular a sum rule which can be derived 
independently of perturbation theory. Specifically, we 
estimated the spin-density structure function for large 
energy transfers to using a quasi-bremsstrahlung ampli- 
tude (Fig. [l]). For small u>, the corresponding neutrino 
scattering rate diverges as 1/uj 2 due to the intermediate 
nucleon going on-shell. Because the true differential scat- 
tering cross section must be finite for all w, and motivated 
by considerations of multiple scattering, we advocated re- 
placing 1/ui 2 by a Lorentzian l/(w 2 + T 2 /4) where T is of 
order r CT , but is adjusted so that the structure function 
obeys the sum rule. 
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FIG. 1. Neutrino-nucleon scattering in the presence of an 
external spin-dependent potential for the nucleons. The po- 
tential can arise from bystander nucleons, a pion condensate, 
the walls in the nuclear bubble phase, or some abstract exter- 
nal force. 

Meanwhile, Sawyer [jl2| has published an explicit treat- 
ment of the cross-section reduction based on more tra- 
ditional perturbative techniques. In addition to the 
quasi-bremsstrahlung graphs of Fig. [I] he includes wave- 
function and vertex renormalizations to elastic scatter- 
ing. The leading correction in the nucleon scattering po- 
tential, V, is the interference between zeroth and second 
order amplitudes, an example of which is shown in Fig. |^. 
These terms diverge, behaving as 5(uj) in the absence of 
nuclear recoil. However, Sawyer points out, the sum of all 
order \V\ 2 contributions yields a total vN cross section 
which is finite, but reduced from the vacuum value. 
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FIG. 2. One diagram representing the interference between 
zeroth and second order scattering amplitudes that leads to a 
wave-function renormalization of the incident nucleon for vN 
scattering. 

Motivated by Sawyer's work, we show how the diver- 
gence of the quasi-bremsstrahlung process represented by 
Fig. [l] can be rigorously controlled by virtue of our exact 
sum rule of the spin-density structure function without 
assuming any specific modification of its form, Lorentzian 
or otherwise, and without calculating the renormalization 
terms explicitly. However, even though our Lorentzian 
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ansatz is not needed to obtain the perturbative cross- 
section reduction effect, it nevertheless yields the correct 
limiting value because this ansatz incorporates the sum 
rule explicitly. Put differently, we may implement the 
sum rule by an explicit ansatz for the low-to behavior of 
the spin-density structure function, or we may use the 
sum rule in an abstract sense. Either way, in the pertur- 
bative limit the final result agrees with the one found by 
Sawyer Q| even though the path of derivation is entirely 
different. Our novel technique has the added benefit that 
after the nature of the perturbative region has been un- 
derstood, the nonperturbative regime may still be stud- 
ied using our proposed Lorentzian modification or some 
other related ansatz. 

In Sect. II we use the structure-function formalism to 
derive the perturbative limit of the average axial-current 
neutrino-nucleon scattering cross section. In Sect. Ill we 
consider nucleons interacting with an external classical 
potential. In this generic example the relationship be- 
tween the perturbative bremsstrahlung matrix clement 
(Fig. and the cross-section reduction becomes particu- 
larly transparent and allows for a direct comparison with 
Sawyer's (jl2| result. 

In a dilute medium where the perturbative approx- 
imation is justified, the most important practical con- 
sequence of nucleon spin fluctuations may not be the 
mild cross-section reduction, but a new mode of energy 
transfer between neutrinos and the nuclear medium p0| . 
This energy exchange is enabled by the nontrivial fre- 
quency dependence of the spin-density structure func- 
tion and thus is specific to spin fluctuations; spin-spin 
correlations do not contribute. Indeed it is plainly vis- 
ible from the bremsstrahlung nature of the underlying 
matrix element (Fig. that neutrinos can transfer en- 
ergy to nucleons above and beyond the standard nucleon 
recoil effect. 

Complementing the numerical expression of Ref. |lC| ] 
we derive in Sect. IV an analytic expression for the av- 
erage energy transfer per collision. This perturbative re- 
sult is relevant for conditions around the neutrino sphere 
in a SN and thus for the formation of neutrino spectra. 
Sect. V is given over to discussion and a summary. 



II. AVERAGE NEUTRINO SCATTERING RATE 

A. Low-Density Limit 

The impact of nucleon spin fluctuations on neutrino 
scattering rates is most easily understood in the long- 
wavelength limit (see Ref. || for a discussion) which has 
been employed in virtually all previous papers dealing 
with neutrino opacities, or neutrino pair and axion emis- 
sivities, in SN cores or old neutron stars. In this limit, 
the momentum transfer between neutrinos and nucleons 
is neglected. The axial-current scattering cross section 
may then be written as 
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where E\ and e 2 are the initial- and final-state neutrino 
energies, Gp is the Fermi constant, and the neutral- 
current axial weak coupling constant in a dilute medium 
is Ca ~ +1.37 and —1.15 for protons and neutrons, re- 
spectively hj. 

For simplicity we focus on an isotropic, nonrelativis- 
tic, nondegenerate medium of baryon density tib, tem- 
perature T, and a single species of nucleons. In this case 
the function S(lo) is the dynamical spin-density structure 
function in the k — > limit laJlT 
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where (<x(t) • <x(0)) is the autocorrelation function for the 
nucleon spin operator cr(t) = J e? 3 x cr(x) at time t. Here 
<t{x) = (x)Ttp(x), ijjix) is the nucleon field (a Pauli 
two-spinor) and r is a vector of Pauli matrices. The 
expectation value (. . .) is taken over a thermal ensemble 
so that detailed balance S(oS) = S(—uj) e"/ T is satisfied. 
We note that our definition of energy transfer is positive 
for energy given to the medium. 

In order to derive an average scattering cross section 
we consider nondegenerate thermal neutrinos which we 
take to follow a Maxwell-Boltzmann distribution; the dif- 
ference to a Fermi-Dirac distribution is inessential for the 
present discussion. Therefore, we consider the quantity 
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With the dimensionless energy transfer x = (si — s 2 )/T 
and after one explicit integration one finds 

r°° dr ~ 

(<ta) = ctt —S(x)(2 + x+±x 2 )e- x . (4) 
Jo 2vr 

Here, a T = f C\G 2 T 2 while S{x) = T S{xT) is the 
dimensionless structure function. In vacuum the nucleon 
spins do not evolve, yielding a constant autocorrelation 
function and thus S{x) — 2-k8{x). Then (a a) = &t where 
/ dxS(x) = g has been used. 

For reasons that will soon become apparent, we con- 
centrate not on a direct calculation of the average cross- 
section at finite density (<ja), but rather on its deviation 
from the vacuum cross section 8 (a a) = (pa) — <Jt or 



8(<t a ) 

(7 J- 



= -1 



dx 
2^ 



S{x) (2 + x- 



\x 2 )e 



(5) 



The crucial step is to express the r.h.s. as a common 
integral over S(x). To this end we use the normal- 
ization J_ °° S(x) dx/2ir — 1 which obtains if the spins 
of different nucleons evolve independently; otherwise an 
additional correlation term would appear on the r.h.s. 
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LO|jm], a possibility to be addressed in Sect. II. B be- 
low. By virtue of detailed balance this sum rule is 



dx 
2^ 



S(x) (l + e- x ) = 1. 



(6) 



Replacing —1 in Eq. (||) by the negative of Eq. (||) yields 
(5(<ta) f°° dx 



2n 



S(x)G(x), 



(7) 



where G(x) = -[(2 + x + f x 2 ) e~ x - (1 + e~ x )] or 



G{x) = l-{l + x+\x 2 )e- x . (8) 

This function is shown in Fig. |^. It expands as G{x) — 
| x 2 + 0(x 3 ) for small x, approaches 1 for large x, and is 

always positive. Because S(x) is also a positive function, 
we find that the average cross section in the medium is 
indeed always suppressed by spin fluctuations. 
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FIG. 3. Function G(x) as denned in Eq. (|). 

Thus far S(uj) has been the nonperturbative but un- 
known structure function. However, what can be cal- 
culated in the framework of perturbation theory is an 
expression Sbremsfw) based on the "bremsstrahlung" or 
"medium excitation" amplitude Fig. |l|. In our previous 
works ]5|,|| , we showed that SWorns^) diverges for small ui 
asw" 2 , a behavior which is generic for all bremsstrahlung 
processes — for the electromagnetic case see Jackson [[TJJ . 
We may then write 



SbremsM = "4 S {^/ T ) 
0J Z 



1 



J ' T for uj < 0, 
for uj > 0, 



(9) 



where s(x) is a nonsingular even function with s(0) = 1. 
The quantity T a , defined as the coefficient of the uj~ 2 
singularity of Sbremsfw), is physically interpreted as the 
spin-fluctuation rateJj 



*One may consider the limit of a classical spin vector s(t) be- 
ing kicked by a random force at a rate r co n . If the spin changes 



It should now be clear why we have calculated the de- 
viation 5{oa) rather than {a a) itself — the x 2 behavior of 
G(x) compensates for the singularity in Sbrcms- We are 
therefore free to substitute Sbrems(^) f° r S(uj) and are as- 
sured of a finite answer for 5((ta}- Further, if we accept 
that in a dilute medium the true S(uj) is well-represented 
by 5'broms(^) for uj > T a then 
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~Z~ S\, TeTns {x) G(x) 
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(10) 



is the desired perturbative result. 

With our representation Eq. (^|) the cross-section re- 
duction is to lowest order in 7 CT = T a /T 



6(<ta) 

CJ'f 



' 2-k 



dx — -~— s(x). 



(11) 



Taking for simplicity the classical limit s(x) = 1 we find 
5(<ta) 



5 TV 

6 2tt' 



(12) 



Once more, these results put in evidence that "f a is the ex- 
pansion parameter which defines the perturbative regime. 

We may estimate the error due to using SWems^) in 
Eq. @ instead of the full S(x) in Eq. @). If the true 
S(x) is given by Sbrems(x) to lowest order in 7cr so that 
§(x) - S hlems (x) = 0{jI) for x > 7<T , then f™[S(x) - 
S hrems (x)]G(x)dx/27r = 0( 1 2 ). 

This implies that the lowest-order cross-section reduc- 
tion effect represented by Eq. ( |ll"| ) will be found by any 
assumed functional form S , ap prox(w) for the true S(oj) if 
SapproxM agrees with S brcms (uj) to 0(y%) for uj > T a . 
Any such function which is normalized can be inserted 
into Eq. (Q) and will then yield Eq. ( 1 1 ) up to an error of 
0(7^). Further, any such function, even if it is not nor- 
malized, will yield this result when inserted into Eq. (Q) 
where the normalization condition has been reshuffled 
into the function G(x). To lowest order in 7^, the cross- 
section reduction effect is independent of the detailed 
structure of the true S(uj) in the neighborhood of uj = 0. 



abruptly by a random amount As in a given collision (which 
is thus assumed to be "hard") and if subsequent spin orienta- 
tions are uncorrelated one finds Shrc-ms(uu) = Y a /{uj 2 + 1^/4) 
with r CT = r co ii ((As) 2 )/(s 2 ) |||||. This justifies identifying 
r CT with an average spin rate of change or a spin-fluctuation 
rate. In the classical limit of hard collisions one has s(x) = 1, 
while for general interaction potentials s(x) is more compli- 
cated. Quantum corrections introduce the detailed-balance 
factor, and cause s(x) to be a decreasing function for large x, 
as discussed for the case of electromagnetic bremsstrahlung 
by Jackson H|. The same conclusion is inferred from the 
f-sum rule for S(ui) 0. 
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B. Spin-Spin Correlations 

A crucial step in the above analysis was use of the 
sum rule in Eq. (^), appropriate for a medium of un- 
correlated nucleons. However, in a real nuclear medium 
the nucleon spin fluctuations are typically caused by a 
spin-dependent interaction among nucleons. Inevitably, 
this will cause correlations between different spins so that 
the r.h.s. of the sum rule Eq. (|J) is 1 + C^a) where in 
a dilute medium |C(7 <7 )| <C 1. It follows that G{x) re- 
ceives an additional contribution — C(j a )(l + e~ x ) and 
8{(ja)Ipt one of order C("f a ). Here 5{(Ta)/&t is to be 
calculated with the full S(x) not St, ienis (x). If C(7cr) is 
of order 7™, then the correction to the cross-section shift 
from considering spin-correlations is also of order 7™. 

We may next use the above estimate of the error in- 
curred by using SbremsOc) rather than the true S(x). 
Then, if C(7 CT ) oc 7", with n > 1 the cross-section de- 
viation calculated from S\> Tenss {x) in Eq. ( |To| ) is to lowest 
order independent of spatial spin-spin correlations. 

For example, if the nucleon-nucleon interaction poten- 
tial is written as in Ref . [|ll| , and if the correlation length 
scales as j a , then we expect C(j a ) to be of order 7^, 
in which case the low-density limit of the cross-section 
change is well described by Eq. (|To|) . 



C. Comparison with the High-Density Behavior 

We next compare the low-density limit thus derived 
with our more general previous expression. In a dense 
medium (j a > 1) the detailed structure of S(uj) for low 
energy transfers matters. In the past we have advocated 
a Lorentzian form 

Sapprox W = S{(J T) X < 

uj< l +T z /A [_1 for lj > 0, 



a plausible intermediate result between the naive and 
lowest-order perturbative results. 



(13) 



where for a given r CT one chooses T such that S, 



approx 

is normalized. This ansatz is motivated by a heuristic 
argument [8| and by the classical limit which obtains for 
oj <C T 15j . Equation ( p^ ) naturally approaches the 
appropriate limit for low densities. 

In Fig. fj we show {oa)/ot for s(x) = 1 as a function 
of 7o-. The dotted line marks the "naive" constant cross 
section which obtains when spin fluctuations are ignored 
entirely. The dashed line represents the perturbative re- 
sult according to Eq. (|l^); for j a > 7.5 it yields complete 
nonsense (a negative scattering cross section). The solid 
line marked "Lorentzian" was obtained with the above 
ansatz for S ap p IOX (ui). The dashed line is its tangent at 
the point 7 CT = so that indeed the Lorentzian ansatz 
yields the same perturbative limit as the direct calcula- 
tion in Sect. II. A where the sum rule was implemented in 
an abstract sense rather than by a specific ansatz for the 
low-w behavior of S(oj). The Lorentzian ansatz yields 




Lj 1 1 I 1 1 1 I 1 1 1 I 1 1 *j 

2 4 6 8 

7a=r„/r 
I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 

.5 1 1.5 

p [10 14 g/cm 3 ] 
(for T=30MeV) 

FIG. 4. Average axial-current neutrino-nucleon scattering 
cross section as a function of the spin-fluctuation rate in the 
classical limit where s(x) = 1. The dotted line refers to a 
naive calculation where nucleon spin fluctuations are ignored 
entirely. The dashed line refers to the low-density perturba- 
tive expansion which agrees with Sawyer's |l2| related result 
if one expresses the spin-fluctuation rate IV in terms of the 
assumed nucleon interaction potential (Sect. III. A). The solid 
line marked "Lorentzian" arises from the ansatz Eq. (|l^ ) while 
the one marked "Top-Hat" arises from Eq. @. The second 
horizontal axis expresses 70- in terms of the physical density 
p for T = 30MeV, using the OPE interaction potential ac- 
cording to Eq. ([l5]). The correspondence between 7 CT and p is 
meaningful only in the perturbative regime where 7 CT < a few. 

The overall shape of the Lorentzian line in Fig. [I| is 
determined by the bremsstrahlung wings of S a p Vmyi {<jj) 
together with the sum rule. In order to test how sensitive 
it is to the assumed \cw-u) shape we have considered a 
second ansatz of the form 



S-. 



approx 



M = r CT x 



for < oj < u>o, 
~ 2 for uiq < lo. 



(14) 



Of course, for oj < we have the detailed-balance factor 
e w / T as in Eq. (|l3|), and in general there is a function 
s(lu/T) which we take to be equal to 1 for the purpose 
of illustration. For a given choice of r CT the frequency luq 
is determined such that S'approx(w) fulfills the sum rule. 
The cross-section reduction derived from this "top-hat" 
ansatz is shown in Fig. ^. It has a common tangent at 
7 CT = with the dashed line and the Lorentzian curve, 
again confirming that in the perturbative limit the de- 
tailed low-w shape of S(oj) does not matter. For large 
7 CT the deviation from the Lorentzian curve is relatively 
small. Therefore, it appears that even in the nonpertur- 
bative regime the cross-section reduction is dominated 
by the bremsstrahlung calculation in conjunction with 
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the sum-rule, and not by fine points of the low-w shape 
ofS(cj). 

In this discussion we have used the classical brems- 
strahlung limit of hard collisions where s(x) = 1. Quan- 
tum corrections alone require that s(x) is a decreasing 
function of x for large x, and the same conclusion is 
reached on the basis of Sigl's f-sum rule H]. Further, the 
detailed large- a; behavior depends on the short-distance 
behavior of the assumed NN interaction potential. It is 
evident from Eq. (|ll|) that the detailed functional form 
of s(x) will determine the slope of the curves in Fig. ^ at 
7cr = 0. However, any change in slope is independent of 
the manner by which S(x) has been adjusted to satisfy 
the sum rule. Thus, the three curves remain tangent to 
each other, although their slope will in general be differ- 
ent from the value — Itt- determined for s(x) = 1. 

To express the spin-fluctuation rate in terms of the 
physical density and temperature we need to assume a 
specific model for the cause of the spin fluctuations. Tak- 
ing nucleon-nucleon interactions modelled by a one-pion 
exchange (OPE) potential as in our previous papers one 
finds for a single species of nucleons 



o-.OPE 



5/2 
l N 



8.6 MeV Pl3 T^\ 



(15) 
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^ 2e 1 2e 2 2E 1 2E 2 V ' V ' 
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x27TS(E 1 +e 1 -E 2 ~e 2 ). (16) 

Here, n c is the number density of classical scattering cen- 
ters, fi is the occupation number of the initial-state nu- 
cleon, k is the momentum absorbed by the external po- 
tential, Pauli blocking factors are ignored for all particles 
because of the assumed nondegeneracy, and the neutrinos 
have been ignored in the momentum S function because 
of the long-wavelength approximation. The expectation 
value (|A4| 2 ) is understood to include the averaging of 
classical ensemble variables on which the external poten- 
tial might depend. The perturbative structure function 
Sbrcms(^) is obtained by comparing Eq. ([l6]) with Eq. (Q). 

In order to derive the matrix clement we use the axial 
part of the weak interaction Hamiltonian and an exter- 
nal classical potential for the nucleon spins. The most 
general form for the potential in Fourier space is ||16|| 



V(k,«r,s) = U {k) + U s (k)a-s 

+ U T {k) (3<r • ksk- er-s) 



(17) 



where a„ = (f2mN/m 7T ) 2 /Air 



pion fine structure constant, pi% = p/10 13 gcm 3 , Ti 



15 with / k, 1.0 is the 

gem" 3 , Tio = 

T/10 MeV, and m N = 940 MeV is the nucleon mass. The 
pion mass has been neglected. Taking T w 10 MeV as a 
typical value for SN conditions, one concludes that the 
dividing line between high and low density is roughly 
given by 10 13 gcm~ 3 or 3% nuclear density. 

We stress that r CT; oPE is in itself a perturbative result 
and thus will be a reasonable representation of the true 
r CT only if T a ope/T ^ 1- Therefore, in the high-density 
regime Eq. ( |l5| ) cannot be used to translate an assumed 
spin-fluctuation rate into a corresponding physical den- 
sity. We have previously argued that the true j a in a 
nuclear medium never exceeds a few fllOlllfl. 



III. MATRIX-ELEMENT REPRESENTATION 



A. Perturbative Cross-Section Reduction 



The perturbative structure function S , b re ms(w) is cal- 
culated from the quasi-bremsstrahlung process shown in 
Fig. |l| so that one may represent the cross-section re- 
duction 5 (a a) directly in terms of its matrix element M. 
The translation is most easily achieved by considering the 
differential scattering cross section. Denoting the four- 
momentum of the in- and outgoing nucleon with (E\, pi) 
and (E 2 ,p 2 ), respectively, we find 



where k = |k|,k = k/fc,sisa classical spin vector of 
length 1 associated with the external potential, and er is 
the nucleon spin operator. Here, Uq is a spin-independent 
potential while Us and Ut represent a spin-dependent 
scalar and tensor force, respectively. After some algebra 
one finds 



E 



(i-Mi 2 ) 



ri2 r<2 



2e 1 2s 1 2E 1 2E 2 2to 2 



\Us(k)\ 



(1 - ±c 12 ) 



-2Re[U T {kM(k)] (cica 



|ci 2 ) 



\U T (k)\ (2-cica 



|c 12 ) 



(18) 



where c,- (i = 1 or 2) is the cosine of the angle between the 
direction of neutrino i relative to k, while ci 2 refers to the 
angle between the two neutrinos. We have averaged over 
the external spin directions s with an assumed isotropic 
distribution. 

The interaction Uo(k) does not contribute because it 
leaves the nucleon spins unchanged. This leaves us with 
the scalar and tensor force Ug(k) and C/t(&)> respectively. 
If the classical scatterers are substituted by the nucle- 
ons themselves only the tensor term survives because the 
scalar term conserves the total spin of two colliding nu- 
cleons and thus does not cause spin fluctuations Q . 

Expression ( |l8| ) reveals explicitly the lu~ 2 divergence 
of Eq. ( |l6| ) which thus cannot be integrated to yield 
a total cross section. However, following the steps of 
Sect. II. A we can derive a convergent expression for the 
medium-induced change of a\ which denotes the total 
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axial-current scattering cross section for a fixed initial- 
state energy e±. In the structure-function language it is 
the ds-i integral of Eq. ([!]) or equivalently 



47T 



^H(ei-^) 2 e(£, 



(19) 



In vacuum <7i 



induced change ba\ = <j\ 
6a! f + °° <Lc 



(3C^G|/47r) e\ so that the medium- 



1 



2tt 



(ei - c^) 2 6(ei - w) 



(20) 



Then we may proceed as before and replace —1 by an 
integral over the structure function by virtue of its nor- 
malization so thatf] 



5<ji 

a 1 vac 



( ei -^) 2 e( £l -^) 



p 2 
"1 



- 1 



(21) 



As before, the integrand varies effectively as S(uj)uj 2 for 
small ui because the term linear in uj switches sign at 
the origin. Therefore, to lowest order we may substitute 
S(ui) — + Sbrcms(^), provided we interpret the remaining 
integral by its principal part. 

£brems( w ) is obtained by comparing Eq. (|l|) with 
Eq. ( |l6| ) and using Eq. (p"8|). After performing the du>, 
dp2, and (IQ2 integrations we arrive at 



2 n c f d 3 Pi d 3 k 



5<j\ 

= 3 ^ J (2tt)3 (2tt) 3 

|2 ,,,|2 



X 



|«7s(fc)r + 2|t7 T (fc)| 



./1 



^) 2 e( £l -w) 



1 , 

(22) 



where the energy transfer is ui — — (2pi • k + k 2 )/2mjv. 



B. Comparison with Sawyer's Result 

As mentioned in the introduction, Sawyer jl2| has dis- 
cussed a cross-section reduction due to the interaction 
of the target nucleons with bystander particles. He does 
not provide an immediate physical interpretation of his 
calculation, but we believe that in essence he has studied 
the same effect that is the topic of our paper, namely the 
scattering version of the Landau-Pomeranchuk-Migdal 



^In this form one can easily see that for small si the cross 
section actually increases. For example, si = leads to 
a vanishing vacuum cross section while in the medium it 
is (3C\Gp/4n) J_ dujui 2 S(u>)/2n or by detailed balance 

(3CiG|/47r) / °° dcjuj 2 S(cj)e-" /T /2ir. 



effect. However, his formal approach is quite different 
from ours. 

The optical theorem implies that calculating the total 
neutrino scattering cross section amounts to a calculation 
of the imaginary part of the neutrino forward-scattering 
amplitude /p on nucleons. Sawyer uses analyticity con- 
straints |17J for /o to recognize that the total cross sec- 
tion should be finite order by order in a perturbative 
expansion in powers of the nucleon interaction potential. 
Further, he observes that this result holds even though 
individual contributions to /o have infrared singularities 
from on-shell intermediate states. One type of V 2 contri- 
bution to fo comes from interference between the zcroth 
and second order scattering amplitudes. Fig. ^ shows 
such a contribution which may be interpreted as a wave- 
function renormalization of the incoming nucleon. Sim- 
ilar terms would renormalize the outgoing nucleon wave 
function, or provide a vertex correction. The other type 
of V 2 contribution to Im(/o) is given by a phase-space in- 
tegral over the square of the amplitudes shown in Fig. [l]. 
These terms correspond to the quasi-bremsstrahlung in- 
elastic scattering process, and also diverge as discussed 
above. Sawyer's main point is that the divergences in 
these two types of terms must sum to a finite result. 

Sawyer [pjj has worked out several examples which il- 
lustrate this approach. Specifically, the cross-section re- 
duction represented by his Eq. (10) is very similar to our 
Eq. (|H). However, our Eq. ( fj2|) has not yet been aver- 
aged over initial state neutrino energies, Sawyer has used 
bystander nucleons to provide the potential so that his 
expression for the energy transfer takes account of the by- 
stander recoil, and he has studied vn — > pe~ scattering 
rather than vN — > Nv so that the proton-neutron mass 
difference appears. Further, he has used a scalar poten- 
tial which is explicitly isospin dependent so that the role 
of our spin fluctuations is played by isospin fluctuations 
in his case. 

The main point of agreement is the structure of the 
term in square brackets in Eq. (^2|). Both Sawyer's 
Eq. (10) and our Eq. ( p2| ) diverge if one considers the first 
term in square brackets independently from the — 1. In 
our derivation, the —1 effectively represents the sum rule 
of the nonperturbative S(u>) for which we have substi- 
tuted 5brcms(w) after the two terms have been combined. 
In Sawyer's approach, the —1 corresponds to the wave- 
function renormalization of the elastic scattering rate. 
Our interpretation of the agreement between these re- 
sults is as follows. 

In effect, Sawyer has calculated the perturbative ap- 
proximation (uj) to second order in V. Recall that the 
zeroth-order approximation is S^(uj) — 2u5{u)) because 
Sawyer also uses the long-wavelength limit where nucleon 
recoil effects are ignored. In this limit, any nonvanishing 
power of S {2) (uj) at uj ^ must arise from the quasi- 
bremsstrahlung amplitudes of Fig. [I] so that inevitably 
S^ 2 '(uj) = 5b re ms(w) for w ^ 0. Sawyer's renormalization 
terms modify only the elastic channel uj = so that his 
complete result amounts to S^(oj) — Sbrcms(^) — A5[iS) 
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where A is an infinite integral expression. Of course, 
S^ 2 '(u)) is highly singular and thus unphysical at u = 
in the sense that in the neighborhood of w = it does 
not provide a representation of the differential scattering 
cross section Eq. ([!]). However, (lu) is legitimate as an 
integral kernel to calculate the total cross section Eq. 
The agreement between Sawyer's and our results shows 
that the second-order perturbative calculation yields an 
expression for A such that S^(u>) fulfills our sum rule. 

In essence, then, Sawyer's calculation amounts to 
showing explicitly that the renormalization terms not 
only cancel the low-w divergence of St,i ema (u>) , but in- 
deed cancel it in such a way that SW(uj) fulfills the sum 
rule. The renormalization terms are an explicit second- 
order manifestation of the information embodied in our 
sum rule. In our derivation, we have shown the sum rule 
to be a general nonperturbative property of S(u>). There- 
fore, once we have calculated SWems^) we can handle its 
Iow-oj divergence either by an abstract application of the 
sum rule, or by an explicit ansatz for the physical behav- 
ior of the true S(lo) near lo = 0. Either way we do not 
need to calculate the renormalziation terms explicitly. 

Although Sawyer's and our approaches are equivalent 
in the low-density limit, they are not equivalent when one 
considers the high-density case. There, a perturbative 
expansion makes no sense as higher-order terms exceed 
the lower-order ones. However, by making use of the 
sum rule, and exploiting the physical insight that S(uj) 
should have a width of order r CT JyJ and possess a hard 
bremsstrahlung tail for u> 3> r^, we have the basis for a 
reasonable model of the high-density regime. 

To summarize, our derivation is based on representing 
interaction rates by virtue of current correlators which 
allow for a direct transition to the classical limit. There- 
fore, our approach allows for an intuitive interpretation of 
the cross-section reduction as a temporal spin-averaging 
effect. Moreover, because we know on general grounds 
that the sum rule Eq. (g) must be fulfilled, we do not 
need to worry about a calculation of the various infinite 
second-order corrections to the elastic scattering rate. In 
our derivation, the only required Feynman-graph evalu- 
ation is that of the "medium excitation term" of Fig. [I]. 
Finally, our derivation allows for a clear and physical 
identification of the dimcnsionlcss parameter 7^ which 
defines the perturbative expansion. Sawyer's technique, 
on the other hand, represents a more familiar method- 
ology if one approaches the problem as a perturbative 
expansion in powers of V. 



IV. ENERGY TRANSFER 



As stressed in Ref. [Q, the most important effect of 
nucleon spin fluctuations may be that they allow for a 
new mode of energy transfer by the quasi-bremsstrahlung 
process shown in Fig. [j]. The relevant figure of merit is 
the average energy transfer per collision (Ae)brcms ° r 



/ °° deieje El/T » J °° de 2 e\{e 2 - e 1 )S{e 1 - e 2 ) 
J °° de 1 e\e-^/ T - J °° d£ 2 e\S{E 1 - e 2 ) 



(23) 



Here, T„ is the temperature of the neutrinos which are as- 
sumed to follow a Maxwell-Boltzmann distribution while 
the nucleons are characterized by T. In Ref. jl7J this 
expression was evaluated numerically on the basis of the 
Lorentzian ansatz for the structure function. 

However, in the dilute-medium limit one can also de- 
rive an explicit expression. We first note that in the 
numerator and denominator one can each perform one 
integration explicitly so that 



(Ae) brcmg _ / °° dx S(x) x F {x) (e~* - e~^) 
r °° dx S(x) Fp(x) (e~ x + e-P*) 



T 



(24) 



with F fJ (x) = 1 + \fix + j2 P 2 * 2 and P = T / T v 

In the dilute-medium limit we may use to lowest order 
Six) — > Sbrems^) in the numerator, while in the de- 
nominator S(x) 2tt 5{x) because the medium-induced 
change of the cross section is itself of order 7 CT . With the 
representation Eq. (|^) we find 

(Ae) brcms = T a s(x) ^ (e - _ e -/J«). ( 25 ) 

Jo 27T X 

For the classical limit of hard collisions where s(x) = 1 
this is 



(Ae) 



-7+6/3 + /3 2 + I21n/3 
24^ ' 



(26) 



This is to be compared with the average energy transfer 
by nucleon recoils, (Ae) recoi] = 30 /T 2 T 2 /m N g§. 

Therefore, the ratio between the two is 



(Ae) 



brcms 



r a m N 



fl 2 



-7 + 6/3 + /5 2 + 121n/3 



(Ae) 

recoil 

T 2 r 720tt(/?-1) 

_ r a TOAr / 1 , 7(0-1) 



T 2 V 36tt 



1447T 



(27) 



Therefore, the importance of the "inelastic" mode of en- 
ergy transfer exceeds that of recoils if j a > 36nT/mN- 

We note that the quasi-bremsstrahlung process of 
Fig. [I] has a standard counterpart where neutrino pairs 
are absorbed or emitted. We define a rate of energy trans- 
fer in this channel, normalized to the average neutrino 
scattering rate in analogy to the above discussion. By 
virtue of Ref. the result can be expressed like Eq. ( |25| ) 
with Fp(x) = /3 s a; 5 /1440. The efficiency of energy trans- 
fer relative to recoil effects is 



(Ae)pa 



r*m N 2 (/3 2 + l)(/3 + l) 



T 2 



3600tt 



(28) 



(Ae) rcco ii 

Therefore, the quasi-bremsstrahlung process of Fig. [j] is 
approximately a factor of 25 more important than pair 
processes. 
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V. DISCUSSION AND SUMMARY 



Bartol, support by DOE grant DE-AC02-78ER05007 is 
acknowledged. 



To summarize, we have studied the neutrino-nucleon 
scattering cross section taking into account nucleon spin 
fluctuations. The effect of random spin fluctuations is 
to reduce the cross section in a manner similar to the 
LPM reduction of photon bremsstrahlung by multiple- 
scattering effects. We have derived perturbative results 
in terms of a lowest-order spin-density structure func- 
tion, and also in terms of the squared matrix element of 
neutrino-nucleon scattering in the presence of bystander 
nucleons or more general external spin-dependent poten- 
tials. In this form, our result agrees with Sawyer's Jl^j 
related finding. The low-density limit is unique unless 
there are unexpectedly strong spin-spin correlations. 

While we have focussed on neutral-current processes, 
similar conclusions would obtain for charged-current col- 
lisions as stressed in Refs. [5 12 . 



The explicit low-density results are theoretically inter- 
esting, but their practical significance is limited. It is 
obvious from Fig. [| that a plausible extrapolation into 
the high-density regime vastly differs from the pertur- 
bative result for j a = T a /T > 1 which implies that the 
perturbative result cannot be trusted for densities greater 
than a few percent nuclear. While T^ope overestimates 
the true IV at nuclear density, in a SN core one has val- 
ues for the true 7 CT of order a few, perhaps as large as 10 
0,0. Therefore, the neutrino opacities in the inner 
SN core cannot be treated by perturbation theory alone. 
Near the neutrino sphere, corresponding to 7^ = 0(1), a 
perturbative treatment is roughly justified, but the cross- 
section reduction is small (a few percent) and thus not 
overly significant. 

Near the neutrino sphere, the most important practi- 
cal consequence of nucleon spin fluctuations is likely to 
be the inelastic or quasi-bremsstrahlung mode of energy 
transfer. With Eq. @ and taking T = 5MeV as a 
typical neutrino-sphere temperature it is found to com- 
pete with standard recoils for 7 CT > 0.5. As this value is 
representative for conditions around the neutrino sphere, 
we confirm that the inelastic mode of energy transfer is 
about as efficient as recoils and thus may be important 
for the formation of neutrino spectra JTo| ] . 
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